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Abstract. The spherical-box approach is extended to calculate the resonance 
parameters and the real part of the wave function for single particle resonances in 
a potential containing the long-range Coulomb interaction. A model potential is taken 
to demonstrate the ability and accuracy of this approach. The calculated resonance 
parameters are compared with available results from other methods. It is shown that 
in the presence of the Coulomb interaction, the spherical-box approach works well for 
not so broad resonances. In particular, for very narrow resonances, the present method 
gives resonance parameters in a very high precision. 
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1. INTRODUCTION 

The investigation of continuum and resonant states is an important subject in quantum 
physics. For the theoretical determination of resonant parameters (the energy and 
the width), in addition to many approaches which are based on the connection of 
the unbound states (including resonances and antibound or virtual states) and the 
poles of the S matrix, several bound-state-like methods have been developed, e.g., 
the complex scaling method (CSM) [U El El SI El El El El E] , the complex absorbing 
potential (CAP) method [TQl [Tl] , the analytical continuation in the coupling constant 
(ACCC) approach P2 Q31 Q31 E51 G51 QH QS QSl ED] , and the real stabilization method 

(RSM) [2HE21E31EIlE51EElEZlEBlESllin]]. 

In the real stabilization method, the Schrodinger equation (or the Dirac equation 
in relativistic models) of the system in question is solved in a basis [21] or a box [21] 
of finite size, thus a bound state problem is always imposed. The RSM is based on 
the fact that the energy of a "resonant" state is "stable" against the change of the 
size of the basis or the box. Many efforts have been made in order to calculate more 
efficiently resonance parameters with the RSM pi E3 EHl EZ1 EE1 EH]. One of the 
effective ways to do so is the so called "spherical-box approach" [21] in which the energy 
and width of a resonance in finite-range potentials are determined from the variation of 
the discrete positive energies with the radius of the box. The spherical-box approach for 
short range potential has been vigorously proved [31]. In a recent work [30], the single 
neutron resonances in atomic nuclei were investigated by combining the RSM and the 
relativistic mean field model [321 [33] . 

In many atomic, molecular and nuclear processes, the long-range Coulomb 
interactions play important roles. Therefore it is necessary to develop the spherical- 
box approach for resonances in the presence of the long-range interaction and check its 
validity and accuracy. In the present work, we extend the spherical-box approach in 
order to include the influence of the Coulomb force. 

The paper is organized as follows. In section El we give the formalism of the 
spherical-box approach for resonances in the presence of the Coulomb interaction. The 
numerical details, the results for a model potential and discussions are given in section [3j 
Finally in section H] we summarize our work. 

2. Formalism 

We deal with a central field problem. The radial Schrodinger equation reads (in atomic 
units) 



1(1 + 1) 



2r 2 



+ V Q (r) + - tt fc (r) =EV k (r), 



(1) 



where Vo(r) is a finite-range potential and E = k 2 /2. Under the box boundary condition 
the continuum is discretized and one is left with a bound state problem. When the size 
of the box 1Z is large enough, the energies of bound states, if there are any, do not change 
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with 1Z. In the continuum region, there are some states stable against the change of 
the size of the box, i.e., the energy of each of such states is almost constant with the 
variation of 1Z; such stable states correspond to resonances. It should be noted that 
the spherical-box approach is essentially different from the "variable phase approach" 
where the potential in question is amputated of its part extending beyond a distance f 
and one needs to solve a series of scattering problems for different f [34J . 

In the spherical-box approach, the resonance energy is determined by the stability 
condition, 

d 2 E 



^ = 0, (2) 

E E>y 

the corresponding box size is labeled as 71, i.e., = E(7Z). 

The width can be evaluated from the stability behavior of the positive energy state 
against the change of the box size around TZ [21] . When r is large enough, the finite- 
range potential Vb(r) vanishes and \I/fc(r) satisfies, 

1 In 

*fc(r) oc — sin(/cr - — — 7 ln(2£x) + r]i), (3) 
KT a 

where rji is the phase shift. — / ~f\n(2kr) comes from the long range Coulomb interaction 
with 7 = Z/k. The box boundary condition ^k(7Z) = gives 
In 

kTZ- — --i \n(2kK) +r}i = nix. (4) 
Thus the derivative of the phase shift with respect to the box size reads, 

dVl L 7 W/ ^+^^+ln(2^)^. (5) 



:)(' 



dK V V2ETZJ V \/2EdTZ J ' dTl 

Around an isolated resonance, the energy E and the phase shift rii(E) satisfy the 
following relation, 

r/2 



Vi(E) — Vi,pot(E) + tan 1 ( 



E Ey 



(6) 



Under the assumption that the phase shift from the potential scattering t]i )POt (E) varies 
slowly with respect to the box size, i.e, dr]i^ ot /cVZ ~ 0, one derives the formula for the 
width, 



r -v 7 ^ 
2 ~ 



(7) 



Note that the above formula goes back to equation (8) in [24] when the Coulomb 
interaction is absent, i.e., Z = 0. 



3. A model problem 



In order to check the ability and the accuracy of the spherical-box approach for 
resonances in potentials containing a long-range interaction, we solve a model problem 
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and compare our results with those predicted by other methods. We choose the following 
potential, 

V{r) = V r 2 e- r + -, (8) 

which has been extensively investigated for Z = and/or Z = —1 [351 I2H [26], [361 13"T] . 
We take Vo = 7.5 in the present work. 

3.1. Numerical procedure 

The Schrodinger equation (DO) with the potential (ED is solved under the box boundary 
condition by using the shooting method [381 Chapter 18] with the forth order Runge- 
Kutta algorithm. For Z = 0, the energy calculated from the shooting method converges 
with the step size decreasing and the relative error reaches to within 10 -8 with a step 
size 5r = 0.001. For Z = —1, a step size 5r = 0.0001 gives the relative accuracy < 10~ 8 
in the energy. In the following we shall present results for the potential ([H]) with Z = —1. 
The step size Sr = 0.0001 will be used in order to get a high accuracy comparable to 
the results given in the literature. 

3.2. Resonance parameters 
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Figure 1. The energies of s wave states calculated in a spherical-box of different sizes. 
The inset shows the first avoid crossing which is zoomed in. The red dotted curve 
represents the potential ([8|) around the barrier. 

We first present results for s states. In figure [1] the potential V(r) given in ([8]) and 
the energies of s states versus the box size are presented. By examining the figure, two 
resonances can be found with energies about 1.8 and 4.0 respectively. The former is 
very narrow as is indicated by the sharp avoid crossings. The latter, however, is above 
the barrier and very broad. 

The resonance parameters evaluated at different stable regions are given in table [1] 
for the very narrow resonance. For comparison, the results obtained from the complex 
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Table 1. The energy and width of the first s wave resonance evaluated at different 
avoid crossings from Q and ([7]). All quantities are in atomic units. 



1Z 




r 


5.0163 


1.780 


531 


212 


6.302 


9xl0~ 


5 


7.0547 


1.780 


525 


482 


8.764 


7x10- 


-5 


8.7509 


1.780 


524 


837 


9.344 


1x10- 


-5 


10.3918 


1.780 


524 


706 


9.507 


7x10- 


-5 


12.0226 


1.780 


524 


661 


9.553 


6x10- 


-5 


13.6542 


1.780 


524 


606 


9.566 


2x10- 


-5 


15.2878 


1.780 


524 


620 


9.569 


6x10- 


-5 


16.9243 


1.780 


524 


629 


9.570 


7x10- 


-5 


18.5633 


1.780 


524 


635 


9.571 


1x10- 


-5 


20.2045 


1.780 


524 


634 


9.571 


3x10- 


-5 


m\ 


1.780 524 536 


9.571 9x10- 


-5 


m 


1.780 5 




9.58xl0~ 5 





methods are also included [361 EZ]. With TZ increasing, the resonance energy converges 
to a stable value very fast. For this narrow resonance, the spherical-box approach gives 
very precise energy and width which are comparable to the exact values given in [37J. 
One achieves a seven-significant-digit accuracy for the energy at the third avoid crossing 
with TZ = 8.7509 and a four-significant-digit accuracy for the width at the ninth avoid 
crossing with TZ = 18.5633. These high precisions are very encouraging. One can then 
safely use this spherical-box approach in studying narrow resonances with computational 
efforts much less than complex calculations. 

Table 2. The energy and width of the second s wave resonance evaluated at different 
avoid crossings from ([2]) and ([7]). All quantities are in atomic units. 



TZ 


E-y 


r 


3.4784 


4.101 765 300 


0.582 256 450 


4.9943 


4.053 255 228 


0.804 538 485 


6.2534 


4.021 498 436 


0.939 289 481 


7.4361 


3.995 348 518 


1.024 872 144 


8.5939 


3.970 907 320 


1.083 354 259 


9.7483 


3.945 799 270 


1.101 418 144 


10.9117 


3.917 591 013 


1.144 053 591 


12.0986 


3.881 431 163 


1.198 720 256 


13.3808 


3.807 212 130 


1.329 057 297 


m\ 


4.101 494 946 


1.157 254 428 



The second s wave resonance is quite broad. In table [2] we list the resonance 
parameters obtained for this resonant state at different stable points TZ. When 
TZ > 13.3808, one can not find a stable region satisfying the condition d 2 E/d7Z 2 = 
in the present numerical accuracy. At the last observed avoid crossing, the resonance 
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energy deviates from the exact value by 8% and the width by 15% [37J . This may 
be used as a guide to estimate the accuracy of the spherical-box approach for broad 
resonances. 

A third s wave resonance was predicted to be lying at 4.66 by the complex 
method [37]. It is very broad with a width 5.34. There is no hint for this resonance 
from figure [TJ This implies that such a broad resonance is beyond the ability of the 
spherical-box approach. We note that for very broad resonances, the assumption that 
the phase shift from the potential scattering t]i^ ot (E) in ([6]) varies slowly with respect 
to the box size may not hold. This could be one of the reasons why the present method 
does not work well for broad resonances. 
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Figure 2. The energies of p wave states calculated in a spherical-box of different 
sizes. The red dotted curve represents the potential ([5]) together with the centrifugal 
potential. 

The energies of p states are plotted as a function of the box size TZ in figure [2J The 
effective potential V e s(r) = V(r) + 1(1 + l)/2r 2 with V(r) given in ([8]) is also shown 
in the same figure. Only one resonance is found which is close to the barrier and the 
energy is around 3.8. This state is broader than the first but narrower than the second 
s wave resonances. The parameters for this resonance are given in table [31 Both the 
energy and the width converge well with TZ increasing. At the sixth stable point with 
TZ = 9.6945, the relative deviation of the energy from the converged value (obtained 
from the last avoid crossing at around 19.8229) is less than 0.1% and that of the width 
less than 1%. 

Figure [3] shows the E ~ TZ plot for the d states. There is almost no pocket in 
the effective potential. But one can still find slightly stable regions in the energy range 
4.5~5.0 in the first several E ~ TZ curves. Indeed we could find stable points fulfilling 
the condition d 2 E/d7Z 2 = at the first six avoid crossings for the d state. Thus the 
resonance parameters could be calculated and they are listed in table [U At TZ = 9.3164, 
the relative deviation of the resonance energy is within 3% and that of the width is 
about 8% compared with the values obtained at the previous stable point TZ = 8.1438. 
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Table 3. The energy and width of the p wave resonance evaluated at different avoid 
crossings from ^ and f7j). All quantities arc in atomic units. 



n 




r 


3 


3342 


3.8665 





1401 


4 


9740 


3.8501 





2101 


ti 


2430 


3.8471 





2388 


r 


4205 


3.8457 





2519 


8 


5641 


3.8450 





2581 


9 


6945 


3.8445 





2610 


10 


8201 


3.8442 





2608 


ff 


9441 


3.8440 





2616 


13 


0681 


3.8437 





2620 


f4 


1925 


3.8431 





2623 


15 


3174 


3.8433 





2625 


1G 


4430 


3.8431 





2626 


17 


5691 


3.8429 





2628 


18 


6958 


3.8427 





2629 


19 


8229 


3.8425 





2631 
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Figure 3. The energies of d wave states calculated in a spherical-box of different 
sizes. The red dotted curve represents the potential © together with the centrifugal 
potential. 

Similar to the case for the second s wave resonance, one can not find further stable 
behavior in the region 1Z > 9.3164. 

No higher partial wave resonances are found in the present study. 

3. 3. Wave functions of the resonances 

Since the schrodinger equation ([!]) is solved within a spherical-box, we can only get 
the real part of the wave function for each resonant states. The wave functions for the 
resonances found in the present work are shown in figure HI The behavior of the wave 
function is consistent with the width of the resonance. For the first s wave resonance 
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Table 4. The energy and width of the d wave resonance evaluated at different avoid 
crossings from ^ and IJTj). All quantities arc in atomic units. 



K 




r 


3.2624 


4.8689 


0.7135 


4.6913 


4.8063 


0.9744 


5.8924 


4.7518 


1.1486 


7.0239 


4.6992 


1.2779 


8.1438 


4.6391 


1.3920 


9.3164 


4.5409 


1.5120 



0.8 
0.6 

s °- 4 

^ 0.2 
0.0 
-0.2 
4.0 
3.0 
§ 2.0 
1.0 
0.0 
-1.0 



Figure 4. The radial wave function 'J'(r) of the first and the second s wave resonances 
(lower panel) and of the p and d wave resonances (upper panel). Note that when 
r > 5.0, the wave functions for the two s resonances are multiplied by 10. 

which is quite narrow, the wave function is almost completely localized inside the range 
of the finite range potential Vo(r) = 7.5 r 2 e~ r . But for the broad s resonance, the radial 
wave function inside the potential barrier is much depressed and it oscillates very much 
in the asymptotic region. For the p wave resonance, one still finds a localization feature 
in its radial wave function, but not so prominent as the case of the first s resonant state. 
The wave function of the d resonance is similar to that of the second s state (note that 
the scales in the upper and lower panels are different), except that the former vanishes 
at the origin. 

The real stabilization method is based on the fact that a resonance is more or 
less localized. Therefore the energy of a resonant state is only weakly affected by the 
variation of the size of the spherical-box. This has been shown in [24"1 [30] as well as in 
the present work. Next we show the convergence of the real radial wave function with 
respect to the box size. The radial wave function of the p wave resonance is shown 
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Figure 5. The radial wave function of the p resonance within a box of different size. 
Note that the wave function is multiplied by the radial coordinate r in order to see 
clearly the asymptotic behavior. 



in figure [5j In order to see clearly the asymptotic behavior, the wave function \&(r) 
is multiplied by the radial coordinate r. With TZ increasing, more nodes appear and 
the wave function inside the potential barrier decreases only slightly. This reveals that 
the wave function of a narrow resonance is affected little by the box size 1Z due to the 
localization property. 

4. Summary 

The spherical box-approach, which is one of the effective implementations of the real 
stabilization methods for single particle resonances, is extended to the case in which 
a long-range force such as the Coulomb interaction plays a role. The formalism is 
presented and the numerical realization is fulfilled for this approach. 

We take a model potential ([H]) as an example to demonstrate the ability and the 
precision of this approach. It is shown that in the presence of the Coulomb interaction, 
the spherical-box approach still works well for narrow resonances. In particular the 
present method can give resonance parameters in a quite high precision for very narrow 
resonances. The energy and width also converge reasonably fast with the box size 
increasing. Within the spherical-box, the real wave functions of these resonances are 
obtained. The localization behavior of the radial wave function is consistent with the 
width of a resonance. The convergence of the radial wave function with the increases 
of the box size is studied, which shows that the wave function of a narrow resonance is 
also "stable" against the changes of the box size. 

Acknowledgments 

We would like to thank Professor Shun- Jin Wang for helpful discussions. This work was 
partly supported by the National Natural Science Foundation (10705014, 10775004, 



Spherical-box approach for resonances in presence of Coulomb interaction 



10 



and 10875157), the Major State Basic Research Development Program of China 
(2007CB8 15000) and the Knowledge Innovation Project of CAS (K JCX3-SYW-N02) . 
The computation of this work was supported by Supercomputing Center, CNIC, CAS. 

References 

[1] Kato K 2006 J. Phys: Conf. Ser. 49 73-78 

[2] Kruppa A T, Lovas R G and Gyarmati B 1988 Phys. Rev. C 37 383-389 

[3] Rcinhardt W P 1982 Annu. Rev. Phys. Chem. 33 223-255 

[4] Ho Y K 1983 Phys. Rep. 99 1-68 

[5] Moiseyev N 1998 Phys. Rep. 302 212-293 

[6] Gyarmati B and Kruppa A T 1986 Phys. Rev. C 34 95-102 

[7] Kruppa A T, Heenen P H, Flocard H and Liotta R J 1997 Phys. Rev. Lett. 79 2217-2220 

[8] Moiseyev N and Cederbaum L S 2005 Phys. Rev. A 72 033605-8 

[9] Arai K 2006 Phys. Rev. C 74 064311-7 
[10] Riss U V and Meyer H D 1993 J. Phys. B: At. Mol. Opt. Phys. 26 4503-4535 ISSN 0953-4075 
[11] Mueller I B, Santra R and Cederbaum L S 2003 Int. J. Quantum Chem. 94 75-92 
[12] Kukulin V I and Krasnopol'sky V M 1977 J. Phys. A: Math. Gen. 10 L33-L37 
[13] Kukulin V I, Krasnopol'sky V M and Miselkhi M 1979 Sov. J. Nucl. Phys. 29 421 
[14] Kukulin V I, Krasnopol'sky V M and Horacek J 1989 Theory of Resonances: Principles and 

Applications (Kluwer Academic, Dordrecht) 
[15] Tanaka N, Suzuki Y and Varga K 1997 Phys. Rev. C 56 562-565 
[16] Tanaka N, Suzuki Y, Varga K and Lovas R G 1999 Phys. Rev. C 59 1391-1399 
[17] Yang S C, Meng J and Zhou S G 2001 Chin. Phys. Lett. 18 196-198 
[18] Cattapan G and Maglione E 2000 Phys. Rev. C 61 067301-4 

[19] Zhang S S, Meng J, Zhou S G and Hillhouse G C 2004 Phys. Rev. C 70 034308-8 
[20] Guo J Y, Wang R D and Fang X Z 2005 Phys. Rev. C 72 054319-8 
[21] Hazi A U and Taylor H S 1970 Phys. Rev. A 1 1109-1120 
[22] Fels M F and Hazi A U 1971 Phys. Rev. A 4 662-674 
[23] Fels M F and Hazi A U 1972 Phys. Rev. A 5 1236-1249 

[24] Maier C H, Cederbaum L S and Domckc W 1980 J. Phys. B: At. Mol. Phys. 13 L119-L124 
[25] Taylor H S and Hazi A U 1976 Phys. Rev. A 14 2071-2074 

[26] Mandelshtam V A, Ravuri T R and Taylor H S 1993 Phys. Rev. Lett. 70 1932-1935 

[27] Mandelshtam V A, Taylor H S, Ryaboy V and Moiseyev N 1994 Phys. Rev. A 50 2764-2766 

[28] Kruppa A T and Arai K 1999 Phys. Rev. A 59 3556-3561 

[29] Adamson S, Kharlampidi D and Dementiev A 2008 J. Chem. Phys. 128 024101 

[30] Zhang L, Zhou S G, Meng J and Zhao E G 2008 Phys. Rev. C 77 014312-6 

[31] Hagcdorn G A and Meller B 2000 J. Math. Phys. 41 103-117 

[32] Vretcnar D, Afanasjev A, Lalazissis G and Ring P 2005 Phys. Rep. 409 101-259 

[33] Meng J, Toki H, Zhou S, Zhang S, Long W and Gcng L 2006 Prog. Part. Nucl. Phys. 57 470-563 

[34] Calogero F 1967 Variable Phase Approach to Potential Scattering (New York: Academic Press) 

[35] Isaacson A D, McCurdy C W and Miller W H 1978 Chem. Phys. 34 311-317 

[36] Yamani H A and Abdclmoncm M S 1995 J. Phys. A: Math. Gen. 28 2709 

[37] Sorianos S A and Rakityansky S A 1997 J. Phys. A: Math. Gen. 30 3725 

[38] Press W H, Teukolsky S A and Vctterling William T F B P 2007 Numerical Recipes: the Art of 
Scientific Computing 3rd ed (Cambridge University Press) 



